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The rotations of microscopic magnetic particles, magnetosomes, embedded 
into the cytoskeleton are considered. A great number of magnetosomes are 
shown to possess two stable equilibrium positions, between which there occur 
transitions under the influence of thermal disturbances. The random rotations 
attain the value of order of radian. The rate of the transitions and the prob- 
ability of magnetosomes to stay in the different states depend on magnetic 
fleld direction with respect to an averaged magnetosome's orientation. This 
effect explains the ability of migrant birds to faultless orientation in long-term 
passages in the absence of the direct visibility of optical reference points. The 
sensitivity to deviation from an 'ideal' orientation is estimated to be 1-2 de- 
grees. Possible participation of magnetosomes in biological effects caused by 
microwave electromagnetic flelds is discussed. 
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Many birds of passage and other migrant animals overcome annually thou- 
sands of miles and accurately find the places of their seasonal habitats. This fact 
is not completely understood yet. A number of hypothesis have been suggested 
to account for this phenomenon j^l. In particular, the optical orientation and 
navigation based on 'maps' of the terrestrial surface and starry sky are studied. 
Also studied is the navigation along the lines of force of the geomagnetic field, 
which are known to be rigidly bound with the geophysical coordinates of the 
Earth. The latter hypothesis finds the experimental validation in the numerous 
facts of the ability of some biological species among microorganisms, insects, 
fishes, birds, and mammals to orient in a magnetic field (MF) or react to the 
change in its direction with respect to other acceptable reference marks. 

At the same time, there is no recognized explanation for this phenomenon 
yet. Magnetic orientation is a part of the more general problem of the biological 
eflticacy of weak, less than 1 G, magnetic fields (MFs). A brief review of the 
theoretical works in this area may be found in PP, and the detailed discussion 
in El- 

The problem is in that magnetic energy of biologically active molecules in 
the geomagnetic field is very small. It does not exceed the energy of the electron 
magnetic moment in the Earth magnetic field 0.5 x 10~^° erg. This is more than 
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seven order less than that of thermal fluctuations, i.e. kT « 4.1 x 10~^ erg at 
physiological temperatures. It is not clear, how could such a small 'signal' cause 
a biological reaction on the thermal 'noise' background. 

However, there are submicron particles, which have magnetic moments. 
They have been found in many living objects, particularly in those displaying 
magnetic navigation ability. They consist of the magnetite mainly. Magnetic 
moment of the particles exceeds the elementary one in 7-9 orders. The energy 
of their turn in a weak magnetic field H is essentially larger than that of thermal 
fluctuations. For single-domain magnetite particles of radius r ~ 10^'^ cm or 
100 nm in the geomagnetic field the energy ^lH « vJH equals approximately 
24kT, where n is the magnetic moment of the particle, v and J « 480 G are 
the volume and the saturation magnetization. 

The cytoplasm near cell membranes is such that the turning of a micropar- 
ticle may serve as a stimulus to cell division or ignite a nerve impulse. For 
example, MF produced by a particle with the afore-cited parameters is of value 
up to 0.1 T near the particle itself and strongly depends on its orientation. 
Therefore, the turn of particle may appreciably change the rate of some chemi- 
cal reactions with the participation of free-radical pairs. 

Particularly interesting are the magnetite particles found in the brain of 
many animals and in human brain. The nerve tissue of the brain is separated 
from the circulatory system by the blood-brain barrier which is impermeable 
for most chemicals. In turn, the circulatory system is separated from the diges- 
tive system. Therefore, relatively large ferro- or ferrimagnetic particles cannot 
penetrate into brain tissue as a pollutant. They are found to have a biogenic 
origin, i.e. they appear over time as a direct result of the crystallization in 
brain matter. Biogenic magnetite particles are often called 'magnetosomes'; 
they were first discovered in bacteria that displayed magnetotaxis [S]. It was re- 
cently shown that magnetic nanoparticles may be produced in DNA complexes 
0]. Biogenic magnetite undoubtedly plays role in the navigation of migrant 
birds ini, insects 13, and in other cases. 

For the explanation of the magnetic navigation, the dynamics of magneto- 
somes was modelled by using the equation of free rotations in a viscous liquid. 
This regimen was assumed to be most favorable from the viewpoint of the mag- 
nitude of possible effects. However, the assumption about free rotations does not 
match the data of electronic microscopy. For example, in some microorganisms, 
magnetosomes are assembled into firm univariate chains, where rotations are 
impossible. The idea was used as well that rigidly bound magnetosome brings 
pressure on a closely set receptor. However, in this case, the energy is trans- 
ferred to the number of molecules simultaneously, so that just a small amount 
of energy relates to each one as compared to the kT. 

It turns out that taking into account the not-too-strict elasticity of the 
medium enables one to describe qualitatively new stochastic rotational dynam- 
ics of magnetosomes that is useful to explain the magnetic navigation. 

This article considers the dynamics of a magnetite particle embedded in the 
cytoskeleton. The latter consists of a 3D net of protein fibers of 6 to 25 nm in 
diameter that include actin filaments, intermediate filaments, and microtubules. 
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Figure 1: Relative arrangements of the vectors of magnetic fields and the 
magnetic moment of a magnetosome. 

The ends of these fibers may be fastened to the membrane surface and to various 
cell organelles. We assume the fibers may also be fastened to a magnetosome 
surface normally covered with a bilayer lipid membrane |^. This fixes the 
position of the magnetosome and constrains its rotation to some extent. The 
stationary orientation of the magnetosome generally does not follow the constant 
MF direction. The balance of the elastic and 'magnetic' torques determines the 
orientation now. The torque m affecting a particle of the magnetic moment 
in an MF H equals m = /i x H. 

Here, putting aside the 3D character of the magnetosome rotations, we con- 
sider the magnetosome's motion in the plane of two vectors: the unit vector n of 
the X-axis, with which the vector of magnetosome's magnetic moment coincides 
in the absence of the MF (equilibrium position, = 0), and the MF vector H, 
Fig.IH 

The Langevin equation for rotational oscillations of the particle is as follows: 



lip + ^Lp + k(f ^ -^H{t) sin{ip — ipo) + ^'{t) , luq — 




where ip is the angular displacement, / is the moment of the particle's inertia, 
7 is the dissipation coefficient, k is the factor of mechanical elasticity resulting 
from the cytoskeleton fibers' bending, £^'{t) is a stochastic torque with the cor- 
relation function {t)^' {t + At)) = 2jKTS{At), while ujq is the eigenfrequency, 
and ipo is the MF direction. Then, we assume the quantity of fibers fastening 
the magnetosome to the cytoskeleton may vary from particle to particle and 
a significant number of magnetosomes are mobile enough to markedly change 
their orientation in the geomagnetic field. This means the mechanical elasticity 
due to the fibers' bending is of the same order as or less than the magnetic 
elasticity k < « 24kT. For magnetite Fe304 particles with the substance 
density p « 5.2 g/cw? and radius r « 10"^ cm, we derive a value loq in the 
order of 10^ rad/s. A resonance, however, is not possible since the inertia forces 
are much less than viscous forces: Iujq <C 7. Hereafter, the inertia term in the 
equation of motion may be ignored. 

The idea of this work is to study the dynamics of a magnetosome fixed into 
a visco-elastic cytoskeleton and predominantly oriented in a direction opposite 
to that of a constant MF. We assume further ipQ = n — -q. The angle 77 should 
be read as the off-course angle, i.e. the azimuth deviation from an 'ideal' or 
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Figure 2: The potential function of a magnetosome at the off-course an- 
gle 77 = and 77 7^ 0. Uq is the potential barrier height, Ui is 
the change of the barrier height with an off-course deviation; the 
elasticity parameter equals a = 0.8. 



reference direction, which is determined by the mechanical bonds that fasten 
the magnetosome in an averaged position with regard, for example, to animal's 
cranium. In other words, in the frame of reference of the geomagnetic field, 
r] is the deviation of animal's orientation from an 'ideal' one in parallel to the 
geomagnetic field vector. 

For small angles 77 the equation of motion takes the form: 'yip + kip = 
^Hsm{ip) + iJLHrjCos{ip) + [t) . With the designations 

k ulH „ 2kT , ^ 

^-llH^ '^-'^ ^^Hh 

the equation is reduced to 

^ + d^U{ip,T)^VDaT) (3) 



with the potential 



U 77) = cos((p) + -ip^ -11 sin((p) . (4) 



Here ^(t) is the centered Gaussian process of unit variance (the identity 5{at) = 
6{t)/\a\ is used). 

The potential energy of a magnetosome in terms of ^iH is shown in Fig. |2| 
As is seen, for not too large angles at a < 1 there are two stable equilibrium 
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positions tp± and the unstable one y'o = 0. Due to thermal disturbances, there 
are transitions from well to well, the random turns being of significant value 
around 2 rad. 

Consider the joint influence on a magnetosome of a random torque ^(i) and a 
magnetic signal related with the deviation of the organism's orientation from the 
reference one. The magnetic signal, which varies directly with ry, results in the 
change of the potential function, Fig. [21 The state of the magnetosome oriented 
in the direction of the absolute minimum of the potential function becomes a 
preferred one. The ratio of the probabilities of the magnetosome to locate in 
the states ip± equals 

p- (JU\ 



^=exp^2-J, (5) 

where 5U — 2Ui is the potential difference of the equilibrium points. This 
implies: 

The 'signal' proportional to the difference of the probability and its equilibrium 
value 1/2 equals 

1 Ui 

s = ~ — . 

^ 2 D 

Since the 'noise' is the equilibrium value 1/2, the signal-to-noise ratio in this 
case is equal to 

Rsn = 2^ . (6) 

The quantities Uq, Ui of the potential Q have no exact analytical presentation. 
Here we derive them as the expansions over the parameter 1 — a, which is 
assumed to be a small one: 

ipl^6{l^a), C/o = ^(l-a)2, ^ 6T^^{l-a) , C/"(0) = a-1, C/"(^±) = 2(l-a) . 

The minimally detectable angle of deviation from the reference course follows 
the equation 

-Rsn = 1 . (7) 

From this equation, substituting the derived value Ui in ©, we arrive to the 
formula: 

"^-'^ = wW=T) • 

This quantity is shown in Fig. O as the function of the elasticity parameter a. 
As is seen, the maximum sensitivity takes place at small values a, i.e. for 'softly' 
fastened magnetosomes. However, arbitrary small values of a make no physical 
sense. 

The expression and then (jSJ are valid for the equilibrium probability 
distribution. This means the changes in the potential have to occur more slowly 
than the relaxation to a statistical equilibrium. One can separate relaxations 
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Figure 3: The minimally detectable off-course angle depend upon the elas- 
ticity parameter. 

within each potential well and between wells. For small values a, when the 
potential barrier is high, the relaxation time is determined mainly by the well- 
to-well transitions. The character time here is the mean first passage time (the 
Kramers time, see for example 



where Tor is a character time of the reorientation of animal, the periods of 'hunt- 
ing' around the reference course. Assuming tor > 1 s and taking 7 « Aniyr^ « 
20 X 10^^^ erg-s for the damping coefficient of the rotations of magnetosome 
with the radius 10^^ cm in a liquid with viscosity v « 10~^ g/cm-s (water), 
the condition is fulfilled if a > 0.65. Therefore, from Fig. O is seen that the 
sensitivity to the course deviations is about 0.03 rad or 1.7 degree. 

We will note, if the signal-to-noise ratio in {Tj) was n times greater, the 
minimally detectable angle would be n times less. 

To explain the observable sensitivity of animals to very weak geomagnetic 
variations, in ni| there was suggested that ensembles rather than single magne- 
toreceptors detect MF changes. Given that, the averaging of the signals coming 
from magnetoreceptors, which cause the signal-to-noise to enhance, occurs in 
the cerebrum. With averaging, the signals of regular rotations are summed; for 
chaotic rotations, what is summed are just the squares of the chaotic signals. 
Therefore, the enhancement of the signal-to-noise ratio is reached proportion- 
ate to the square root of the number of magnetosomes contributing to endpoint 




The equilibrium distribution takes place on condition that 



= tor > T, 

7 



6 



response. We remind that the density of magnetosomes in human brain tissue 
was measured around 5 x 10^, and in brain meninges more than 10^ crystals per 
gram '.J^ ■ Some birds have the concentration of the magnetite crystals in orders 
higher. 

There is another factor of the signal-to-noise enhancement that apparently 
was not reported in literature earlier. If we assume each single bird makes flight 
corrective turns during passage in accord with the flight path of the flock, i.e. of 
the most of birds, then the presence of an effective averaging mechanism follows 
that is taken out of the bird's organism. Evidently, this mechanism reduces the 
level of fluctuations in flight direction in proportion to the square root of the 
number of birds in flock. 

Even without an additional enhancement of the signal-to-noise ratio, the 
retention of the orientation accurate to 1-2 degrees during passage could provide 
its exactitude since the additional course correction may occur due to the optical 
navigation in time periods of good visibility. 

In conclusion, we note the following. (1) The stochastic dynamics of mag- 
netosomes has peculiar features in ac MF (a stochastic resonance of magneto- 
somes) and in the slowly varying geomagnetic and 'zero' magnetic fields as well. 
These features enable one to explain the biological effects of weak extremely- low- 
frequency MPs, geomagnetic storms, and a 'magnetic vacuum'. (2) The Debye 
relaxation in water has a wide maximum in the microwave range of the electro- 
magnetic radiation. The viscosity dispersion of water vito) is related with this 
maximum. High-frequency oscillations of magnetosomes in the microwave field 
result in a change of the cytoplasm viscosity in the layer nearby magnetosome's 
surface. This affects the mean frequency of random rotations of magnetosomes 
and may cause a biological response. 
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